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^ ■ Abstract: We consider the two-brane Randall-Sundrum (RS) model with bulk gauge fields. 

' We carefully match the bulk theory to a 4D low-energy effective Lagrangian. In addition to the 
o I four-fermion operators induced by KK exchange we find that large negative 5* and T parameters are 
o • induced in the effective theory. This is a tree-level effect and is a consequence of the shapes of the W 
^ and Z wave functions in the bulk. Such effects are generic in extra dimensional theories where the 
standard model (SM) gauge bosons have non-uniform wave functions along the extra dimension. 
' The corrections to precision electroweak observables in the RS model are mostly dominated by 
^ ■ 5". We fit the parameters of the RS model to the experimental data and find somewhat stronger 
■ ' bounds than previously obtained; however, the standard model bound on the Higgs mass from 
precision measurements can only be slightly relaxed in this theory. 



1. Introduction 



Theories with extra dimensions might explain some of the outstanding problems of particle physics 
[1-3]. In particular some of these models could shed light on why gravity is so much weaker than 
the other three forces. One of the prominent proposals of this sort is the Randall-Sundrum (RS) 
model [2,3], where the strong warping of the extra dimensions introduces an exponential hierarchy 
between the Planck and the weak scales. There are several variants of this model, depending 
on whether the extra dimension is finite (RSI) or infinite (RS2), and whether or not the gauge 
fields are in the bulk. Each of these models can be interesting for slightly different motivations. 
Here we will concentrate on the case where the extra dimension is finite (so that it solves the 
hierarchy problem), and where the gauge fields are in the bulk. This model could possibly yield 
unification of gauge couplings [4], and also may have a simple physical origin [5,6] via the AdS/CFT 
correspondence [7]. The holographic dual of this theory should be a broken conformal field theory, 
which becomes strongly interacting at low energies and spontaneously breaks the weakly gauged 
SU{2) X U{1) electroweak symmetries. Thus this holographic dual of the RS model with the 
gauge bosons in the bulk is in essence a technicolor-like theory, where the broken CFT replaces the 
technicolor group [5] , and the KK modes of the gauge fields and gravitons would be interpreted as 
bound states of the CFT resulting in the technimesons, analogously to the glueball states appearing 
in the case of ordinary AdS/CFT [8]. In QCD-like technicolor theories the new strong interactions 
introduced to solve the hierarchy problem always generate large contributions to the electroweak 
precision observables [9], in particular there are large contributions to the S'-parameter. However 
the phenomenological importance of (non-QCD-like) approximate conformal symmetry has long 
been emphasized in the technicolor literature [10], where the slowly running gauge coupling is 
refered to as "walking". The difficulty of estimating the value of S in these walking theories is 
also well known [11] since it involves non-perturbative, non-supersymmetric gauge dynamics near 
a non-trivial fixed point. Therefore it is interesting to find out whether there is a non- vanishing 5"- 
parameter in the RS model since it provides us with the first approximately conformal ( "walking" ) 
model of electroweak symmetry breaking where such a calculation can be performed. However, in 
the 5D gravity theory (that is the RS model) the S parameter should not be the effect of quantum 
loops, but rather a purely tree-level effect. The purpose of this paper is to carefully match the RS 
model to an effective 4D description and find the value of S in the effective Lagrangian describing 
electroweak physics in this model. Indeed we find that the wave functions of the W and Z bosons 
are distorted due to the Higgs expectation values on the TeV brane, resulting in different wave 
function and mass renormalizations of the W and Z. The physical consequence of this effect is non- 
vanishing S and T parameters, which we calculate. Our method of finding the low-energy effective 
4D theory is general, and we expect that similar effects will appear in any extra dimensional theory 
where the SM gauge bosons have non-uniform wavefunctions. In addition to these parameters the 
well-known effect of the four-fermion operators generated by the exchange of Kaluza-Klein gauge 
bosons has to be included. The coefficient of these four-fermion operators has been called V 
in [12, 13]. We use a global fit to the most recent precision electroweak data to place a bound on 
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the size of the extra dimension in the RS model, and find bounds that are somewhat stronger than 
those previously obtained. 

The paper is organized as follows: in Section 2 we review the results on gauge propagators 
and wave functions in the RS model that will be necessary to calculate the effective Lagrangian. 
In Section 3, we match the higher dimensional theory to an effective 4D Lagrangian, and evaluate 
the S and T parameters. In Section 4 we first calculate the V parameter and then use these results 
for constraining the parameters of the RS model via a global fit to the electroweak precision 
measurements. We conclude in Section 5, while Appendix A contains the detailed expressions of 
the electroweak observables in terms of S, T and V and the SM input and experimental values 
used for our fit. 



2. The Gauge Propagator and Wavefunctions in the RS Model 



In this section we review the results on gauge propagators and wave functions in the RS model 
[4,12-17] that will be necessary for us to calculate the effective low-energy theory. 



The 5D metric of the RS model can be written in the form 
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^ (_^^2 ^ ^^^^^^ ^^u^ ^ ^2.1) 

for i? < 2; < R'. Here R represents the radius of curvature of the AdS space. There is a Planck 
brane a,t z — R and a TeV brane Bit z — R' which cutoff the space with Z2 orbifold boundary 
conditions. 

The 5D action for the bulk SU{2) x U{1) gauge bosons is given by 
J Jr I 4:9i 4:95 

+}^^Jl_SlG^^{Wl,Wh + WlWl + {W!, - - Sp))l , (2.2) 

where the 5-function mass terms arise from a localized Higgs expectation value (H) — v/2. Since 
5-functions on boundaries require special care, we will take the definition of this term to be the 
limit of having the Higgs localized at a point which approaches the brane. This amounts to a 
factor of two difference in the definition of v'^ from taking the Higgs directly on the brane, but 
has the advantage of making comparison with the calculation via the AdS/CFT correspondence 
simpler [18]. 
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We define the weak mixing angle through s, which represents the bare value of the sin 9w- 



s — 



95 



c = 



95 



(2.3) 



We can diagonalize the action by performing a field redefinition: 

W'^^c^Z^ + A^, B^^-s^Z^ + A^. (2.4) 

The reason behind this unusual form for the field rcdcfintions is that none of the fields W^, B, Z 
or A are canonically normalized, but it is equivalent to the standard redefinition in the canonical 
basis. In our new basis we obtain the Lagrangian: 



'5D 



R' ^ 

d'^x I dz — 

R Z 



2 ' MN ' 

5 

2 



/2 



4 \9i 9'5 



FmnF 



MN 



- ^ ,,, ZMiv^^^ + ^<5(^ - I^)-W^W-^ + ^5(^ - R')-ZmZ 



M 



(2.5) 



In the i?^ gauge where FF5 — Z^ — A^ — 0, the propagator for the bulk gauge boson is given 

by [4], 



(2.6) 



where 



(2.7) 



and AvF satisfies 



^5^|5(z - z') = (9,^ - + q'- \v'gl5{z - R')^)Awiq, z, z') 



with boundary conditions 



1 R 



(2.8) 



(2.9) 



also propagates in a generic R^ gauge, but the coupling to fermions is pseudoscalar, and 
vanishes at zero fermion mass. In addition, because is fixed to be odd under the Z2 (consider 
the Z2 behavior of IVs^), it vanishes on the TeV brane and hence does not couple to matter on the 
TeV brane. For our purposes we will therefore be able to neglect W^. 

For R < [z, z') < R! the Green's function can be constructed by patching together the solutions 
of the corresponding homogeneous equation with z < z' and z > z\ which we refer to as A^< and 
/S.w> respectively: 

Aw = 0{z- z')Aw> + 9{z' - z)Aw< (2.10) 
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Plugging the patched solution into Eq. (2.8) for z' ^ z' ^ R' yields: 



^W<\z=z' — ^W>\z=z' 

d,{Aw> - Aw<)\z=z' = gl^ (2.11) 

Setting z' = R' in Eq. (2.11) and combining with Eq. (2.9) yields the IR [z = R') boundary 
condition for the propagator with a source on the TeV brane: 

A I 2 "^2 2 A I /o 1 o\ 

(yzAw<\z=z'=R' — ~^95 ~ '^Qb'" ■^Aw<\z=z'=R' (2-12) 

The solution has the form 

Awiq, R', z) = Aw<\z'=R' = ^ {awMl^) + PwYiiqz)) (2.13) 

where the coefficients are given by 

AglR'YojqR) -AglR'JojqR) 

- D(q) ' m ^^-^^^ 

and the denominator is 

D{q) = Jo{qR) {AqR'YoiqR') + gyRY^{qR')) - Y^{qR) {AqR' UqR!) + g^R^qR')) . (2.15) 

Setting z = R' gives the propagator on the TeV brane. The coefficients aw and I3w have the same 
denominators, and the roots of these denominators determine the 4D poles of the propagator (when 
the numerators do not vanish concurrently). The n*^ pole corresponds to the n*'* W eigenmode 
, with mass m'^\ We will label the lowest mode by n = 0. It is this lowest mode that 

we would hke to identify with the observed W gauge boson, so we will write M^-* = Mw For 
vR -C 1 this is an "almost zero mode" and we can find the pole analytically by expanding the 
Bessel functions in qR < qR' -C 1. To leading order in the coupling the pole is at: 

The bulk Z propagator is obtained from Aw by taking gl ^ g'l + g'i, which gives: 

Rlog{R'/R)4R'^ ^^-^^^ 



Similarly we can find the wavefunction, ij^w^ eigenmode. The wavefunctions are 

given by: 



,(n). . ^ z MMi-h)Yo{Ml;;^R) - Y,{Mi-^z)MMl-^R) 

^^^^ R'Mm^^r')Yo{mPr)-y^{mI-^r')Mm1^^r) 
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where we have normahzed the wavef unction by tp^^^R') — 1. 

For vR <^ 1 the n = mode is an "almost zero mode" and we can find a simple expression for 
the wavefunction by expanding the Bessel functions in M^R <C 1, M\yR' <C 1. To order we 
obtain: 



A/r2 

i/j^^\z) ^l + ^[z^-R'^- 2zHog{z/R) + 2R'Hog{R'/R)] . (2.19) 



The wavefunction for the Z can be obtained by Mw Mz, while the photon remains massless 
and its wavefunction is simply i/:!^\z) — 1. 



3. The 4D Effective Lagrangian of the RS Model 



Given the wavefunctions we can find the 4D effective action by integrating the 5D action over z. 
We want to match the RS calculation onto an effective 4D theory After integrating out the Higgs, 
the most general Lagrangian for the electroweak gauge bosons (with operators of dimension 4 or 
less) can be written as [9, 19,20]: 



1 1 1 sr 



where = 1 - U'^^, Zw = I - ^ww^ Zz = I - U'zz, ^'-yz, ^wwiO), and 112^(0) incorporate 
the effects of new (oblique) physics beyond the standard model. Although we are only doing a 
tree-level matching calculation we have adopted the standard notation for vacuum polarizations 
to represent the wavefunction renormalizations that arise from classical 5D physics. With our 
conventions / ~ 246 GeV. 

Using the 5D wavefunctions (to second order in masses) we can easily calculate the coefficients 
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of the kinetic and mass terms: 



1^ 



JR ^ 5*5 + 5*5 



9' + 9" 91 + 9'^Jr 



? + (^n..(o) ^^^'^^ [ = + n33(o) (3.2) 



where 



n',: = — (2«- >og(«7«) - 2/?^ + g^) 



"33(0) = J' , 2R" - J—--- + = -^n;, . (3.3) 



64i?'4 \^ \og{R'/R) \og{R'/Ry 

64i?'4 log(i?Vi?) ^ \og{R'/Rf 

Here we have used the leading order results for and Mz, Eqs. (2.16). The corrections to the 
wave function renormalization H']^^ and Ilgg arise from integrating the W and Z wave functions, 
while the contributions to the mass renormalizations Hn and 1133 appear from the 5D kinetic terms 
of W and where a z derivative acts on the wave functions. 

Even though there are no l/(167r^) loop suppression factors, for vR <^ 1, Hn and IIss can be 
treated as small perturbations of the leading terms. Thus a simple convention is to identify the 4D 
bare gauge couplings with the leading terms and the Z^ — l with the subleading 11' terms. This is a 
convenient choice because we want to separate out the new 5D physics from radiative corrections 
by loops of standard model particles. It also ensures that there are no additional corrections to the 
couplings of the W and Z to quarks and leptons. In other words, with this convention the mixing 
angles determined by diagonalizing the 5D action are identical to the bare 4D mixing angles: 

- c^^^. (3.4) 



It is these mixing angles that appear in the quark and lepton gauge couplings. Other conventions 
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are possible, but physical observables are independent of convention. Thus in our convention 

1 /I 1 



+ ^)Rlog{R'/R) 

95 / 



\ = \Rlog{R'/R) 
9 9i 

-^Am = ^^^iog(^V^)- (3.5) 

9^ + 9'^ 95 + 95 

Using this identification of the 4D couplings we then obtain for the other parameters of the effective 
Lagrangian: 



R 



2 



f = -^v' , IlwwiO) = 5'nn(0) , Uzz{0) = {9' + ^")n33(0) . (3.6) 



Note that since the photon is massless it receives no 5D renormalization, so 11^^ = ^^^qq — 0. 
Furthermore since this is a tree-level calculation no new Z — 7 mixing can be induced so H^q — 0. 
We can now use the standard definitions [9,21] for the oblique parameters: 

S = 167r(n^3 - U',q) 

u = i67r(n;i - n^3) (3.7) 



Plugging in our results yields: 

S ^ -47rfR'Hog{R'/R) 
T^-^j'R"logiR'/R), 

[/ = , (3.8) 

where wc have dropped terms which are suppressed by powers of log{R'/R). Note that both 5* 
and T are negative and large (i.e. log enhanced relative to a naive dimensional analysis estimate). 

We can check these results by examining the poles of the propagators directly by expanding 
the denominators for qR <^1. At leading order the W pole is determined by 

At next to leading order we must keep terms that are suppressed by q^R'^ relative to the leading 
terms. The pole is then determined by 

1 i?2 ^RXogiR'/R) ,2^4 
= + ^ 2 + V + Bq^ (3.10) 
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where 



A^^^{2R'Hog{R'/R) + R'-R") , 

B = - A ^ \og{R'/R) + R^- R'"") . (3.11) 



Thus to sub-leading order the pole is at: 

.2 d2„2 ^4 7?2,,4 



M2 ~ ^ ^^""^ f^^/' , i2R!Hog{R!/Rf - 2R:Hog{R!/R) + i?'^) 

^ R\og{R'/R)AR'^ 64i?'Mog(i?Vi?)3 ^ / ; &v / ; ^ y 

- ^ - ^ (2i?" ^og(R'/R? - 2R'^ log(R'/R) + R'^) , (3.12) 

4 64 

which agrees at this order with the effective Lagrangian calculation 

^9'{j + nn(o)) (1 + /n;o ^ + ^n;,) . (3.13) 



4. The Comparison of RS to Data 

In addition to the oblique corrections we have described, once quarks and Icptons arc included in the 
theory with couplings to the bulk gauge bosons they will have additional four-fermion interactions 
beyond those in the standard model due to the exchange of the gauge boson resonances. The eflfect 
of these corrections has been parameterized in ref. [12, 13] by a correction to Gp denoted by V. 
Recall that the contribution to Gi;' from W exchange is [21]: 

To include the effect of resonances we can write the bulk W propagator as a sum over poles: 

Aw(q, R', R') = gl y ^ \ (4.2) 



where A^„ is determined by 



and hence 



7V„= p dz\^^^\z)\\ (4.3) 
Jr 



No^ ZwRlog{R'/R) . (4.4) 
Since we chose V£^(i?') = 1, we then have 



A^(g, R', R') = +9'f: . (4.5) 
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At zero momentum the first term on the right hand side is just — 4-\/2G'f,w, and the remaining 
terms are the additional corrections not coming from the W pole. If we write the correction to 
as Gp — Gf,w{^ + then we have 



V^- L^iq = 0, R', R') + TT^^) + nn(0)) 



4 

where in the last line we have again dropped terms suppressed by powers of \og{Rl /R). 

Thus we see that there are three types of corrections to precision electroweak observables 
in RS models: 5", T, and V. Note that some of these corrections have also been considered in 
Refs. [15,17]. To relate our parameters to observables we use the standard definition of sin^o from 
the Z pole, 

sin^eocos^^o = 4^^, (4.7) 
sin^ do = 0.23105 ± 0.00008 (4.8) 

where [22] q;(M|)-^ = 128.92 ± 0.03 is the running SM fine-structure constant at Mz- We can 
relate this measured value with the bare value in this class of models, 

sin^ ^0 = + ^ (-^^ + aT-v), (4.9) 

which is obtained by considering all corrections to (4.7) in the usual way (see [21]). Also, in the 
RS model we have the simple result that with only the tree-level 5D renormalizations the running 
couplings defined by Kennedy and Lynn [23] which appear in Z-pole asymmetries are the same as 
the bare couplings: 

sl{q') = s\ el{q') = e\ (4.10) 



In addition to the contribution from T, there are further corrections to the low-energy ratio 
of charged- to neutral-current interactions coming from resonance exchange. We will absorb this 
efi^ect into the parameter p^: 

=f + nn(0) fl + V/c^\ _ „ s2 



, .^l + aT+—V . (4.11) 

£+n33(o) V 1+^ ; ^ ^ 

Curiously in the RS model we find that the contributions from T and V cancel, and = 1. 
We will however present the general results for precision observables in Appendix A without 
assuming a relation between T and V, so that our results can be used for more general models. 
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With Eqs. (4.9) and (4.11) it is straightforward to calculate the corrections in a general model 
to precision electroweak observables in terms of S, T, and V. The expressions for the various 
observables together with the SM predictions and experimental results are given in Appendix A. 

The result of a global fit to the 23 observables listed in Table 1 is that for Mniggs = 115 GeV 

R'{\og{R'/R)y/^ < 0.50 TeV"^ (4.12) 

at the 95 % confidence level. Taking \og{R'/R) — 32 (as is often done to naturally explain the 
hierarchy between the Planck and weak scales [2-4]) we have: 

l/i?'>llTeV. (4.13) 

For a value of R' which saturates this bound we have 

S = -0.19 
T = -0.03 

V = 0.00082 (4.14) 

Since the T and V contributions to cancel and for these values the contribution of S to (s^ — 
sin^ ^o) is about 8.6 times larger than that of T and 2.6 times larger than that of V, it is the -S" 
parameter constraint that dominantly determines the bound on R'. 

It is interesting to note that for the Z-pole observables, which do not depend on p*, one can 
absorb the contribution from V into an effective T: 

Teff = r--. (4.15) 

a 

Thus in the RS model, T^s is even more negative than T. One can then use the bounds on S and 
T to estimate the bounds on R', which yields results similar to (4.13). 

In the RS model there is also a hght radion that contributes to precision electroweak observables 
at the loop level (together with loops of the KK gravitons and gauge bosons) . These contributions 
have been calculated separately [24] and are small unless an extra Higgs-radion coupling [25] is 
introduced. With this additional coupling, the radion corrections tend to make S more negative 
and hence only tighten the bound on R'. 

In the SM, the fit to data gets significantly worse when the Higgs mass is raised. The reason is 
that the Higgs contributes positively to S, while the data prefers a small or negative S. However, 
in our case S is negative, so one might think that a larger Higgs mass can be accommodated. 
Unfortunately at the same time the Higgs also contributes negatively to T, as can be seen from 
the approximate expressions [21] 

1 / \ 3 / \ 
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Therefore even though the agreement with S can be improved by raising the Higgs mass, T will 
start to deviate even more. In order to see if the fit can be improved we have repeated it for 
the SM results evaluated at Mh = 300 GeV and 600 GeV. In the case of the SM, for the 23 
observables listed in Table 1 in the Appendix increases by about 11 as increases from 115 to 
300 GeV. If we now turn on the corrections from the RS model, the difference between the minimal 
X^'s for the 300 GeV and 115 GeV Higgs reduces to about 7.6, slightly improving the fit, but still 
outside the 95% confidence region (Ax^ = 6.2) of a two parameter fit (our parameters being the 
Higgs mass and X — pR!'^\og{R' / R)). Hence, the Higgs mass bound is shghtly relaxed but not 
significantly. Assuming a 300 GeV Higgs in turn would relax the limit for R' to 1/R' > 9.0 TeV 
(again assuming log R'/R — 32). For the case of the 600 GeV Higgs the increase in minimum x^ is 
19, and is clearly excluded by a wide margin in the two parameter fit. Assuming a 600 GcV Higgs 
the bound on R' becomes 8.2 TeV < 1/R' <22 TeV. These results are illustrated in Fig. 1. If one 




115 GeV 



f^R'^ log (R'/R) 

Figure 1: The change in in the RS model as a function of R' and R, for three different values of 
Mh = 115,300, and 600 GeV. 



were intent on having a heavy Higgs one could add additional new physics to the model that give 
positive contributions to T [27], but this seems completely ad hoc in the present context. 

The bound (4.13) on R' pushes the W gauge boson resonance masses up to around 27 TeV. 
The graviton KK masses are always heavier than the gauge boson modes [13], for example here 
the lowest possible value is around 46 TeV. Thus in the RS model with bulk gauge bosons there 
will be no possible signatures from resonances in WW scattering to observe at the LHC. The focus 
would have to be on the Higgs-radion sector [26] . 

To the extent that the RS model corresponds to a technicolor like model we find these models 
can be consistent with experiment as long as the resonance masses are dialed up. Unfortunately 
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for technicolor models there was no parameter which changed the resonance mass independently 
of the electroweak breaking scale /. 



5. Conclusions 

We have matched the 5D RS model onto a 4D low-energy effective theory and found large (loga- 
rithmically enhanced) negative contributions to S and T. It is interesting to note that this is the 
first model to naturally produce a large negative value for S [27]; however this is mitigated by the 
fact that the resulting bound on R' from precision electroweak measurements forces a seemingly 
unnatural hierarchy of about a factor of 50 between the Higgs vev and the nominal electroweak 
scale of the RS model 1/E! ^ 11 TeV. This pushes the gauge boson resonances (and the graviton 
KK modes) far beyond the reach of the LHC. 
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Appendix 

A. Predictions for Electroweak Observables 

In this appendix we give the predictions of a general model model with contributions to S*, T, 
and V for the electroweak precision observables. We also give in Table (1) the experimental data 
[22,28] and the SM predictions used for our fit in Section 4. Using the results given in [21,29] as 
well as the low-energy ve couplings: 

Qevij^e ^ ue) ^ 2p^{s'^ - ^) , g^A{ve ^ ve) ^ . (A.l) 



12 



we find the following results: 

Tz = (Tz)sM (1 - 3-8 X lO'^-^ + O.OllT - 1.41/) 

= (/?e)5M (1 ~ 2.9 X 10-^5 + 2.0 X IQ-^T - 0.26V) 

= {R^)g^ (l - 2.9 X IQ-^S + 2.0 X lO'^T - 0.26V) 
R^ = {Rr)sM (1 - 2-9 X 10-^5 + 2.0 X lO-^T - 0.26V) 
ah = {(Th)sM (1 + 2-2 X 10-^5 - 1.6 x IQ-^T + 0.021^) 
Rb = {Rb)sM (1 + 6-6 X 10-^,5 - 4.0 X lO-^T + 0.052^) 

Rc = (Rc)sM (1 - 1-3 X 10-=^^ + 1.0 X io-=^r - o.isy) 

^FB = (^fb)sm - 6-8 X 10-^^,5 + 4.8 x lO^^^T - 0.62^ 
= (^fb)5m - 6-8 X 10-^,5 + 4.8 x lO-^^T - 0.62^ 
^FB = (^fb)sm - 6-8 X lO-^-S + 4.8 x lO^^T - 0.62^ 

Ar(Pr) = (Ar(Pr))sM " 0.0285 + 0.0207 - 2.6V 

Ae(Pr) = (A^(Pr))sM - 0.028,5 + 0.020r - 2.6V 

^FB = {^fb)sm - 0-020'5 + O.OMT - 1.8V 

A'pB = {Apb)sm - 0-016,5 + O.OllT - 1.4^ 

Alr = iALR)sM - 0-028,5 + 0.020T - 2.6V 

Mw = {Mw)sM (1 - 3-6 X 10-^,5 + 5.5 x 10-=^T - 0.71^) 

Mw/Mz = {Mw/Mz)sM (1 - 3-6 x 10~^S + 5.5 x lO^-'^T - O.llV) 

gl{uN vX) = {gliuN J^X)) - 2.7 x 10"^ S + 6.5 x lO^^T - 0.066^ 

gliuN vX) = {gliuN J^X)) + 9.3 x lO'^S - 2.0 x lO'^T + O.IOK 

gevii^e ue) = (gevi^^e ^ z/e))^^ + 7.2 x lO-^S - 5.4 x lO^^T + 0.65V 

geA{iye ue) = {geA^ve -> z/e))^^ - 3.9 x lO'^T - 0.151/ 

Qw{Cs) = {Qw{Cs))gj^ - 0.793,5 - 0.0090T - 95V (A.2) 
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